Topological phases in adiabatic and non-adiabatic driven systems 
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In this work we study the geometrical and topological properties of non-equilibrium quantum 
systems driven by ac fields. We consider two tunnel coupled spin qubits driven by either spatially 
homogeneous or inhomogeneous ac fields. Our analysis is an extension of the classical model in- 
troduced by Berry with the addition of the spatial degree of freedom. We calculate the Berry and 
Aharonov-Anandan geometric phases, and demonstrate the influence of the different field parameters 
in the geometric properties. We also discuss the topological properties associated with the different 
driving regimes, and show that by tuning the different parameters one can induce topological phase 
transitions, even in the non-adiabatic regime. 

PACS numbers: 



Introduction: Geometric phases in quantum physics 
is a fundamental issue which has been addressed in the 
last decades. Since Berry's discovery of geometric phases 
in quantum systems [l| , and their further applications to 
different branches of physics (e.g. condensed matter and 
optics) , geometry, topology and nature seem to be mixed 
at a highly complex level. The characterization of differ- 
ent physical properties in purely geometrical and topo- 
logical terms, such as quantized transport [2] and electric 
polarization|3j, has lead to the understanding of some 
undergoing processes governing the physical features of 
these systems [4]. An important characteristic of geo- 
metrical and topological properties in physical systems, 
arises in their strength against different perturbations. 
Some examples of those are backscattering in edge states 
and fault tolerant quantum computation[5|-|8j. 

Theory of Principal Fiber Bundles (PFB) has also 
played an important role in quantum mechanics due to 
its suitability to describe the underlying geometrical and 
topological properties^ [lj|. Its relation with condensed 
matter was pointed out by Simon [Tl| , and the existence 
of an universal PFB [1 2l . Il3l | , which is a purely geometric 
object, shows that non-adiabatic processes can be used 
for practical purposes in a similar way than their adia- 
batic partners. 

The interplay between geometry and ac-fields is a very 
promising field with a huge activity in the last years. 
The tuning of geometrical properties by the applica- 
tion of ac-fields offers a very interesting scenario, where 
non-equilibrium systems show quantum phase transi- 
tions, and the possibility to create topologically protected 
states[IM3|. 

As Berry demonstrated, the phase acquired during the 
cyclic evolution of a general quantum system can be 
expressed as a combination of two terms, a dynamical 
phase 7d and a geometric phase 7g , being the latter ob- 
tained due to the parallel transportation of the vector 
state through the base manifold [18]. The parallel trans- 
portation is obtained by means of the 1-form connection 
A n := i(n\d\n) (d is the exterior derivative operator, and 



|n) is the instantaneous eigenstate), and the geometric 
phase acquired is given by Jq = J c A n , that only de- 
pends on the path followed in the parameter space. 

These geometric quantities describe the evolution of 
the system as the external parameters are varied, but 
other quantities, such us Chern numbers (topological in- 
variants) can also be defined through the 1-form A n , 
characterizing the topological properties of the whole 
base manifold. 




Figure 1: Schematic figure for the relation between adiabatic 
and non-adiabatic regimes. The function /* maps curves C 
from the Bloch sphere (projective Hilbert CP N_1 ) into curves 
C in the parameter space S 2 = {9,ip}. The U (1) fibers at- 
tached to each point due to the phase invariance of quan- 
tum mechanics are represented by small circles. We consider 
(Qp-W-i w ith N = 2 for the schematic figure (i.e.. CP 1 which 
is the Bloch sphere for a single qubit). 



In the present work we consider an extension of the 
classical model considered by Berry [l|, including the spa- 
tial degree of freedom, inhomogeneous magnetic fields 
and iron adiabatic evolution. The analytical results show 
that the interplay between the spin and spatial degree of 
freedom is remarkably relevant when the ac-field varies 
between different sites. We calculate the geometric phase 
and Chern number for both adiabatic and iron adiabatic 
evolution. This model can be realized in different exper- 
imental setups, and the results can be easily extended to 
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larger size systems such as atomic arrays with a pseudo- 
spin degree of freedom. We demonstrate that adiabatic 
evolution for the ac-field allows tunable geometric phases 
by varying the external field parameters: B (intensity of 
the external field) and <f> (phase difference between sites) . 
We show that the variation of the phase difference <fi, 
which introduces the field anisotropy, strongly modifies 
the Berry phase and the topological invariants. Also that 
non-adiabatic processes increase the tuning possibilities 
with the addition of new topological phases, which only 
appear out of the adiabatic regime. 

Model: The system contains a qubit tunneling be- 
tween two sites (L, R) in presence of and ac-field: 



H(t) = B cos (9) SI +J2 tLRcl 



energies now depend on the ratio between the ac-field 
and the Zeeman splitting through 8 (Figj2]) , and not just 
on the total intensity B = \B\ = \JB^ C + B%. 



(1) 




Y, B sin (9) [cos ( Vi (t)) Si + sin ( w (t)) S* y ] . 

i=L,R 

The coupling between sites is given by t^R, B z is the 
static field which splits the qubit levels, and the ac-field 
B ac (t) couples the different states of the qubit (perpen- 
dicular to B z ). In this Hamiltonian, the field has been 
parametrized according to the angles of a 2-sphere 5 2 = 
{6,<p}, being B z = Bcos{9), B x ac = B sin (9) cos ( Vl (*)) 
and B% c = B sin (9) sin (ipi (t)). This is a very appro- 
priate parameter space for a circularly polarized ac field 
with time dependence ipi (t) = Qt + fa, which is the case 
considered in this work. For simplicity, we consider sym- 
metrical B z and B ac (t) at each site. 

For this Hamiltonian we calculate the geometric and 
topological properties by means of the connection 1-form 
A n . We consider both, the adiabatic and non-adiabatic 
time evolution for phase differences <f> = <\>l~4>r £ {0, 7r}, 
i.e. magnetic fields in phase or in phase opposition re- 
spectively. We show that in the last case, interesting 
topological properties arise. Furthermore, we calculate 
the phase diagram demonstrating that non-adiabatic evo- 
lution leads to the appearance of new topological phases. 

Adiabatic case: The adiabatic case considers slow 
time evolution. This allows to define instantaneous eigen- 
states and neglect the transitions to other energy lev- 
els. The instantaneous energies E n ((f)) can be charac- 
terized by two indexes after the diagonalization n = 
(mi = ±l,m2 = ±1): 



E n (0) 
E n (n) 



B 
~2 

m\B 



(m 1 +m 2 X) (2) 
v/l + A 2 - 2m 2 Acos(0), 



being A = 2tLf>/B. We can see that both instantaneous 
energies are ip independent. The case <fi — reflects the 
classical result obtained by Berry (i.e.. without 9 depen- 
dence), with an extra splitting due to the tunneling tLR- 
However, the case (f> = tt shows a 6 dependence which has 
not been obtained previously. That is, the instantaneous 



Figure 2: Instantaneous energies E n (tt) vs 9 for tLR = 1, 
4> = 7r and B — 2 (dashed lines: E±,± (it) and continuous 
lines: -E T ,± (it) ). 



In the calculation of the instantaneous eigenstates, due 
to the parameter space which is a 2-sphere, we need to 
consider two different charts. This is due to the fact 
that the surface of a sphere cannot be directly mapped 
to a plane. The transformation of the instantaneous 
eigenstates between different charts is given by phase 
factors [10]. This is a signature of non-trivial geometrical 
and topological properties. However, we shall not worry 
about this, because we are interested in the geometric 
phase and the curvature tensor F n = dA n (needed for 
the calculation of the Chern number) , which are globally 
defined for the case of Abelian theories. 

In order to obtain the geometric phase, we calculate, 
using the instantaneous eigenvectors, the connection 1- 
form A" = (A™,Ag). The result shows that the phase 
difference strongly modifies the geometric phase (Fig(3]), 
which is tLR dependent for the case of <f> — tt: 

7£(4> = 0) = jr(l-micos(0)) (3) 
mi (Am 2 - cos (6)) + f m2 (A, 



7S = tt) 



fm 2 (X,0) 



(4) 



being f m2 (A, 9) = y/l + A 2 - 2m 2 Acos (0). Note that 
for 4> = (spatially homogeneous field) Jq does not de- 
pend on iiij. This fact, implies that it remains invariant 
for several tunnel coupled two level systems as long as 
4> = (i.e. the field is spatially homogeneous). Once we 
add the spatial anisotropy through the phase difference 
0^0, Berry's result is modified. It adds an interplay 
between tunneling and ac-field which has not been pre- 
viously analyzed. 

The next step is the calculation of the curvature 2-form 

F n = dA n = (j_ A n ,> 



qq A™ 1 dtp A J, that will be used 
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Figure 3: Geometric phase in the adiabatic limit 70 vs 
for different A (A = (left), A = 0.6 (center) and A = 1.2 
(right)). Dashed (continuous) lines: <j> = n (<f) — 0). Note 
that as A = 2tLn/B is increased, a gap appears for the case 
<f> = n. The value at which the topological phase transition 
occurs is exactly where the gap opens (A = 1). 



for the calculation of the Chern numbers: 
sin (0) 



T7171 
^0 = 



f: 



mi 



mi 



2 

sin (0) 



1 — ATO2 cos (0) 



(l + A 2 -2Am 2 cos (6)f /2 ' 



(5) 
(G) 



The result shows that for = 7r, the original curvature 
(Eq[5]) is renormalized by a factor which depends on t lr 
and 77i2 (Eq[n]). For A = 1, the curvature 2-form inverts 
its value, and at the limit of large tunneling t^R 3> B 
the curvature tends to zero for all 0, evolving towards a 
flat space. The first Chern number is given by: 



C F n 



mi<d (1 — A) for (f>=n 
mi for cf>=0 



(7) 



being 9 (x) the Heaviside step function. Hence, for <fi — 
7T, the system undergoes a transition from a topological 
phase Z (ci ^ 0) to a trivial phase (ci = 0) at A = 1 
by tuning the ratio between the hopping t^R and the 
intensity of the ac- field \B\. 

All results in this section has been obtained under the 
adiabatic assumption, but this is an approximation for 
ideal slow evolution. Therefore a natural question arises: 
How do the previous results change when we consider 
non- adiabatic evolution?. 

N on- adiabatic case: Aharonov and Anandan (A-A) 
proposed a generalization of the Berry phase to non- 
adiabatic processes, the A-A phase 7a-a [HI- As later 
was pointed out, this generalization corresponds to a 
universal connection A in a universal PFB CP°° or A- 
A bundle It establishes a general relation between 
the classification of U (1) PFB and quantum mechanical 
systems [13]. The relation between Berry's and Anadan's 
result is stablished through a function /* (pullback bun- 
dle), that maps curves C defined in CP°° to curves C in 
the parameter space, i.e. C = /* (C), then /* has all the 



information for the computation of the geometric phase 
(the explicit form of the Hamiltonian determines /*). 

We will show below that for the non-adiabatic case, 
the geometric phase 7a-a and the phase diag ram match 
the previous results in the adiabatic limit [20]. 

Considering the present setup in the non-adiabatic 
regime, we shall demonstrate that it can be exactly solved 
(without the adiabatic approximation) , and calculate the 
geometric phase 7a-a and the topological invariants. 

By means of a unitary transformation to a co-rotating 
frame with the ac-field, we obtain a static Hamiltonian 
which can be diagonalized. The unitary transformation is 
given by U (t) = exp {-iClt (S^ + ) }, and the trans- 
formed Hamiltonian H = U [ HU — iUW reads: 

h = y, ( B cos w - ") s i + E 

+ B sin ^ t cos (&) s * + sin (&) s l] ■ ( 8 ) 

i=L,R 

which is independent of t. Note that the time dependence 
of the transformation renormalizes the energy levels by 
Q. 

The diagonalization of the Hamiltonian leads to the 
energies £ n (<j>) : 



B 



-mi — y/l + [i 2 - 2/it cos (0) - m 2 t LR 



- mi -Jl+ A 2 -2A m2 cos(0) (9) 



being \i = fl/B and A„ l2 = (0+ ^mithR) /B. Note that 
in both cases, the limit — > matches the adiabatic re- 
sult. Note also the similarities between EqJ5] and Eq(5] 
where the frequency renormalizes the different parame- 
ters. 

Interestingly, the energies present degeneracy points 
for both 4> = 0, 7T. Our calculations show that the geomet- 
ric phases are also Abelian at degeneracy points due to 
the structure of the Floquet operator, which is in agree- 
ment with the results of [2l|. 

The solutions to the Schrdinger equation are given by 
the transformation of the eigenvectors obtained by 
diagonalizing EqJS] to the original frame, i.e. (t)) = 
U ((p (i)) Therefore, the calculation of the 1-form A v 
gives: 



A 



ip 



ity {t)\d v W{t))d<p = $\S z $)d<p, 



being the A-A geometric phase 7a-a 
i.e. : 



fc 



7a-a(0) = mm 
7a-a<» = m i 7r 



(fi — cos (9)) 



yjl + /i 2 - 2^iCOS (6) 
(A m2 -cos(g)) 
v/1 + A^ 2 - 2A m2 cos(0) 



2nA v , 

(10) 
(11) 
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and the curvature: 



F n (0) 
F n (tt) 



sin (0) 



1 — /Lt COS (0) 



(l + /i 2 -2/icos (0)) 3/2 



(12) 



mi- 



sin (6>) 



(1 



,cos(0)) 



1 



A 2 

m 2 



2A m2 cos 



(*)) 



3/2 



We can see that the geometric properties for = does 
not depend on tLR nor on mi (always a degeneracy for 
states with different mi is present). By contrary, the 
= 7r case presents different curvature for each state, 
but because of the functions A± cannot be independently 
modified by changing the system parameters B, Q and 
tt,R, the geometric phases and curvatures also cannot. 
The Chern numbers can be calculated through the cur- 
vature, and then, obtain the phase diagrams as a function 
of A m2 or /«: 



C? (0) 

C?(7T) 



mi 



(l + sign(l- M )) 



(13) 



mi9(l-|A m2 |) 

mi0(l - |0 + 2m 2 tifl| /B) 



The results obtained for the geometric phase show cer- 
tain interesting features. For the — case, the 
non-adiabatic phases 7a-a present some corrections in /j, 
(which is precisely the adiabatic parameter) to the adia- 
batic case, but they remain degenerate in pairs, and the 
variation of /i modifies all geometric phases simultane- 
ously. Also we must note that 7^_a (4> — 0) fEq llOp . as 
in the adiabatic case (EqJ3J) , does not depend on m.2 nor 
on tLR, meaning that systems with a larger number of 
coupled spins would have the same geometric phase. 

Regarding the phase diagram, the = case just con- 
tains two phases: The topologically trivial and the Z 
phase. The trivial phase appears when we go out of 
the adiabatic regime by increasing the frequency (/x = 
il/B > 1), while the Z phase governs the whole adia- 
batic regime. 

Concerning the = 7r case in the non-adiabatic regime, 
the geometric phases (Eq lll|) depend on two functions 
A TO2 , and hence on both indexes {m\,m2) of the state 
vector. This allows to modify the geometric and topo- 
logical properties differently for each state. This is an 
important difference with the = case, because now 
the geometric phase for each state behaves differently 
when we vary tLR, and B for = tt. This is reflected 
in the phase diagram, where two new topological phases 
appear: (0,Z) and (Z, 0) (first index labels states with 
m 2 = +1 and second index states with m 2 = —1) in 
which only two of the states have non-vanishing Chern 
number. More interestingly, at arbitrary high frequencies 
fl we show that states with non- vanishing Chern number 
can appear (FigQJ. 

Because the functions A± depend both on the system 
parameters, we plot the phase diagram as a function of 



B and f2 for fixed tLR- Interestingly, the (Z, 0) phase is 
not accessible (FigQJ. This is due to dependence on B. 
tLR and ft of the functions A±. Note that for a differ- 
ent dependence on the parameters, the (Z, 0) phase could 
appear. We also observed that the variation of the pa- 
rameter t^B only modifies the width of the (0, Z) phase, 
such that in the limit tLR — > we obtain the phase dia- 
gram for the = case. 




Figure 4: Topological phases diagram vs (B, Q) for fixed 
tLR = 1 and <f> = tt. The phase diagram shows the physi- 
cally accessible topological regions. Note that the adiabatic 
limit is found by choosing Q = with the trivial and fully 
topological phases described above. In addition, for non adi- 
abatic processes (i.e. SI > 0) a new phase (0, Z) appears, 
in which just half of the states are characterized by a non 
vanishing Chern number. 

The existence of phases with non-vanishing Chern 
number implies a non-trivial bundle, and therfore holon- 
omy elements different from the identity. It means that it 
is impossible to avoid the geometric phases in the study 
of these systems. The geometric phases are local, in dif- 
ference with the Chern numbers, and can be different 
from zero even for regions where c% = 0. Our results can 
be important for experiments involving the measurement 
of geometric phases, which can be performed by using su- 
perconducting qubits and quantum state tomography[22] 
among other setups 23]. 

Summary: Summarizing, we analyzed the geometric 
phases and topology of tunnel coupled systems with an 
on-site pseudo-spin degree of freedom, and driven by ac- 
fields. We have considered the case of an electron spin, 
tunneling between two different sites, which is coupled to 
an ac magnetic field that produces transitions between 
the spin up/down levels. The interplay between the spa- 
tial and spin degree of freedom due to the ac magnetic 
field, as well as its dependence with the external field pa- 
rameters B, f2 and (phase difference between different 
sites) is analyzed from a geometrical point of view. 

Interestingly, we found that a system presenting spatial 
anisotropy (in our case due to the phase difference be- 
tween the external ac driving fields in different sites) leads 
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to a complex and rich behavior compared with Berry's 
classical result for a localized spin in the adiabatic regime. 
Our analysis is also extended to the non-adiabatic regime. 
In this case we find geometrical phases which depend in 
a non trivial way on the spatial field anisotropy. The re- 
sults generalize the well known Aharonov-Anandan phase 
for localized spins, resulting in a novel topological phase 
diagram (Fig|4j) . Our results are the basis for the anal- 
ysis of larger systems with spatial periodicity which will 
be the subject of a future work. 

Although the topological phase transition can be dif- 
ficult to be directly measured in our setup, the mea- 
surement of geometric phases could indirectly show the 
topological phase transition (see Fig|3]). Setups involv- 



ing quantum dots in slanting magnetic fields [24j and 
quantum circuits[22| can be used for this purpose. Also, 
the non-adiabatic analysis shows how the topological in- 
variants change out of the adiabatic regime, leading to 
possible applications in non-adiabatic quantum compu- 
tation. 

We acknowledge MAT 2011-24331 and ITN, grant 
234970 (EU) for financial support. A. Gomez-Leon ac- 
knowledges JAE program. 
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